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l/-\ , The specific heat (C m ) and optical birefringence (An) for the magnetic percola- 

tion threshold system Feo.25Zno.75F2 are analyzed with the aid of Monte Carlo (MC) 
simulations. Both An and the magnetic energy (U m ) are governed by a linear combi- 
nation of near-neighbor spin-spin correlations, which we have determined for An using 
MC simulations modeled closely after the real system. Near a phase transition or 
when only one interaction dominates, the temperature derivative of the birefringence 
[d(An)/dT] is expected to be proportional C m since all relevant correlations neces- 

o 

sarily have the same temperature dependence. Such a proportionality does not hold 
for Feo 2sZno 75F2 at low temperatures, however, indicating that neither condition 
above holds. MC results for this percolation system demonstrate that the shape of the 
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temperature derivative of correlations associated with the frustrating third-nearest- 

O 
O 
O 



neighbor interaction differs from that of the dominant second-nearest-neighbor inter- 
action, accurately explaining the experimentally observed behavior quantitatively. 
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Measuring the linear optical birefringence (An) in anisotropic, antiferromagnetic crystals under- 

o" 

going magnetic phase transitions is a powerful way of determining the magnetic specific heat (C m ) 
critical behavior. It has been shownuu that the temperature derivative of the optical birefringence 
[d(An)/dT] is proportional to C m . In many cases, the birefringence technique has provided the 
most precise experimental determinations of universal critical behavior parameters!] in pure and 
randomly mixed and dilute magnetic Ising systems. For the case of the three-dimensional (d = 3) 
random- field Ising model (RFIM), which applies for a dilute antiferromagnet with an applied field 
along the ordering direction, birefringence measurements yielded evidence of a phase transitionB. 
The advantages of the technique are threefold: the technique is typically much easier to employ 



than traditional heat pulse techniques; the effects of concentration gradients inevitably present in 
mixed and dilute crystals can be greatly reduceda; and the phonon contributions to the specific heat 
are greatly suppressed in the birefringence. Since the transition typically varies with concentration 
in mixed and dilute systems, the critical behavior is often masked by the concentration gradients 
quenched into the system during growth. The laser beam used in the optical technique can be 
aligned perpendicular to the gradient, often reducing the gradient effects by an order of magnitude. 
This has been crucial in the study of random-exchange Ising model (REIM) and RFIM systems in 
d = 2 and 3. The virtual elimination of the phonon background has allowed detailed analysis of 
the specific heat in d = 1 and 2 systems. For d = 2 this has allowed a detailed comparisonQ of the 
magnetic specific heat of the pure system and the Onsager solution to the d — 2 Ising model and, 
for dilute systems, a scaling analysis of the destruction of the phase transition by random fields in 
the d = 2 Ising modelQ. For d = 1, the technique has been used successfully to determine the ex- 
change constants in systems where the phonon background overwhelms the magnetic contributions 
to conventional specific heat dataH. 

That C m and d(An)/dT are proportional in certain cases is a consequence of the fact that both 
An and the magnetic energy (U m ) are governed by a linear combination of near-neighbor spin-spin 
correlations!! This proportionality has been experimentally shown explicitly^ for the d = 3 pure 
Ising system FeF% as well as its site-diluted counterparts&jO Feo.g3Zno.oTF2 and Feo.4eZno.54F2 
by comparing directly the pulsed heat and optical data. This proportionality has been demon- 
strated both in zero applied field (H = 0) and in applied fields. The magnetically dilute crystals of 
Fe\- x Zn x F2 have been the most extensively studiedta realizations of the RFIM, an important class 
of systems with randomness and frustration. Such experimental verification has been crucial in this 
controversial field since some authoraI3 have expressed reservations about the proportionality for 
H > 0, i.e. the RFIM case. 

Despite the experimental verification of the proportionality between C m and d(An)/dT, it is 
important to verify the details of the mechanism for the proportionality near phase transitions in 
a specific case. One opportunity to analyze the detailed relationship between C' m and d(An)/dT 
presents itself in the compound Feo.25Zno.75F2, which, with respect to the dominant exchange 



interaction, J2, is at the magnetic percolation threshold, the concentration below which no ordering 
is possible. Throughout the rest of this work we will refer to the percolation threshold concentration 
as the appropriate one if only Ji were being considered. The presence of other interactions make 
the concept of percolation more complicated and they can play an important role in the ordering 
processes even when they are very small. Interestingly, d(An)/dT data for Feo.25Zno.75F2 exhibit 
a changelij in sign near T — 8 K. While the high temperature behavior of d(An)/dT is accurately 
proportional to C m , the sign reversal in d(An)/dT at low temperature is not reflected, of course, 
by the behavior of C m . We present data obtained from both the optical and pulsed heat techniques 
that show the nonproportionality between C m and d(An)/dT at low T, as shown in Fig. 1. We also 
present Monte Carlo (MC) simulation results that can be used to explain why the proportionality 
breaks down between C m and d(An)/dT at low temperatures. U m and An each depend on the near- 
neighbor spin-spin correlations, but with different weights. We show that whereas U m is dominated 
by the second-nearest-neighbor correlation, An is predominantly governed by the third-nearest- 
neighbor correlation that shows a different temperature dependence below 15 K. The results of 
this study are consistent with theory regarding the proportionality between C m and d(An)/dT 
near a phase transition, where the different spin-spin correlations must necessarily have the same 
temperature dependence. We also demonstrate that an applied field has a predictable effect on both 
U m and An. 

The compound Fe\_ x Zn x F% is ideal for RFIM studies. Pure FeFi is well modeled by the Hamil- 
tonian 

I i<j i 

The 5 = 2 system has the interaction strengths J x = -0.069 K, J 2 = 5.278 K, and J 3 = 0.279 K, 
as determined from neutron inelastic scattering measurementso The large crystal-field anisotropy 
D = 9.29 K persists while the exchange interactions remain largely constant as the magnetic spins 
are dilutedllj in Fe\^ x Zn x F2 1 making this an excellent Ising system for all magnetic concentra- 
tions x. Single crystals can be grown for all x with very small concentration gradientsES and with 
superb structural quality. The magnetic interactions are dominated by the antiferromagnetic second- 
nearest- neighbor super-exchange interaction, J2, between the body-center and body-corner ions. All 



other interactions are negligible except near the percolation threshold concentration, where a small 
frustration due to the third-nearest-neighbor interaction between ions of the same sub-lattice along 
the direction perpendicular to the spin-ordering c-axis, J3, becomes importantO. At low tempera- 
tures, the Heisenberg character of the spin-spin interaction is less important and we might expect 
U m to be fairly represented, for H = 0, by the REIM Hamiltonian, 

H = EE JwjSiSj , (2) 

l i<j> 

where Si — ±2, e, = 1 if site i is occupied and zero otherwise, and Ji is the strength of the I th 

nearest-neighbor interaction. C m is then given by a sum of temperature derivatives of spin-spin 

correlations, 

d{U m ) sr^sr^ T d(< S l S J >;) 
l i<j 
where the Ji are assumed independent of temperature. The temperature dependent birefringence 

(An), the difference between the indices of refraction along the spin ordering axis and perpendicular 

to it, depends only on a linear combination of the same correlation functions. The temperature 

derivative yields 

rf(An) v-v- r d{< SjSj >i) 

I i<j 

In general the values of the coefficients 7j associated with An are not related to the respective values 
Ji. Nevertheless, we will find a proportionality between C m and d(An)/dT in a temperature regime 
where the temperature dependence of the contributing correlations < SiSj >; are all proportional. 
This is necessarily the case in the critical region of a phase transition where the correlation length 
grows larger than any relevant interaction length and all spin-spin correlations necessarily have the 
same temperature dependence. We may also find an excellent proportionality between C' m and 
An in the case where one magnetic interaction dominates in both C m and d(An)/dT. Such is the 
caseEl of Afni 7 ^, where the proportionality holds accurately over a very wide temperature range, 
5 < T < 100 K. A lack of proportionality has been observed in pseudo-low-dimensional systemsaH 
where neither of the above conditions holds. The percolation threshold presents another opportunity 
to observe a breakdown in proportionality between C m and d(An)/dT, since small frustrating inter- 
actions can affect spins with few or no neighbors connected by the dominant exchange interaction. 



We will show that this is indeed the case for the percolation threshold sample Feo. 25-2^0.75-^2 • We 
will also show that the behavior is readily interpreted using the results of MC simulations of this 
Ising antiferromagnet. 

The Feo. 25-^^0.75^2 crystal was cut from a boule growrO at the University of California, Santa 
Barbara. For the adiabatic heat pulse technique, the 2.24 g sample was mounted on a thin sapphire 
plate using GE7031 varnish. A small Stablohm 800 wire heater, connected to a four-wire constant- 
power supply, was wound and varnished onto the sapphire plate. A shielded carbon thermometer 
was attached with varnish to the sample and was connected using a four-wire technique to a cur- 
rent ratio transformer resistance bridge. The sample was suspended inside a sample chamber by 
0.0254 mm Be-Cu wires used to connect to the thermometer and heater. While providing good 
electrical connections, the alloy wires provide only a small heat leak from the sample to the copper 
sample chamber. In this way, the heat generated by the thermometer leaks away from the sample in 
a way controlled by a temperature difference, ST, between the sample chamber and the sample. In 
the absence of a heat pulse, the sample temperature can remain constant if ST is properly controlled. 
This is done by controlling ST with a bridge and controller using the sample thermometer and a 
thermometer located in a cold finger, commercially calibrated for H = 0, upon which the sample 
chamber is mounted. A second thermometer in the cold finger is used with a second bridge to 
determine the absolute temperature. A preliminary calibration is performed to determine ST versus 
T such that the sample temperature does not drift in the absence of a heat pulse. The calibration is 
incorporated into computer control algorithms so that the balance is automatically preserved over 
the entire temperature range during a specific heat experiment. The computer applies a pulse and 
determines the resulting change in temperature, from which the specific heat is calculated. The ther- 
mometry sensitivity is approximately 50 /iK. The pulse energy ranged from approximately 3 /xJ at 
the lowest temperatures to 8 mJ at the highest temperatures where measurements were made. After 
collecting specific heat data for the sample and addenda, the sample is removed from the sample 
chamber and the thermometer is fixed to the sapphire plate. The specific heat of the thermometer, 
the sapphire plate, the varnish and the wires is then measured at H ~ 0. This background specific 
heat is subtracted from all of the specific heat data before further analysis. The phonon contribution 



to the total specific heat is approximated by the Debye model, valid at temperatures small relative 
to the Debye temperature. The Debye temperatureEj for pure non-magnetic ZnF 2 is 250 K and that 
of FeFi is 256 K. We weighted the two pure Debye temperatures by the respective concentrations 
to obtain a Debye temperature of 252 K for Feo.25Zno.75F2- The Debye specific heat is calculated 
with this Debye temperature and is subtracted from the data leaving the magnetic component of 
the specific heat (C m ) which can then be used in comparisons with An and MC results. 

For the birefringence (An) measurements, two faces were polished parallel to the spin-ordering 
c-axis. A linearly polarized laser beam (A = 632.8 nm) impinges normally upon the sample with a 
polarization 45° to the c-axis. The beam traverses a distance of 8.05 mm through the sample. A 
0.5 mm pin-hole in front of the sample minimizes the effects of concentration gradients and vibrations. 
The Senarmont technique is used to measure An with a resolution of 2 x 10~ 9 . Calibrated carbon 
resistance thermometers, used for their low field dependence and high sensitivity, yield a temperature 
stability better than 50 /iK. The sample was zero-field-coolcd (ZFC) to 5 K before slowly raising 
the temperature in approximately 0.1 K steps. About 400 s were required to establish equilibrium 
and to measure An at each step. Reasonable variations in the rates of heating and cooling and 
stabilization times had no observable effect on the data. 

The correlation functions for the first three nearest-neighbor pairs have been calculated using 
MC simulations. The magnetic lattice corresponding to the body-centered-tetragonal Fe.\_ x Zn x F2 
lattice is described as two cubic sub-lattices of size L x L x L each, delineated as one dimensional 
arrays bit coded to accommodate large lattice sizes. All of the results reported here were obtained 
with L = 256, corresponding to more than 3.3 x 10 7 sites magnetically occupied with probability 
x. Beginning at high temperature each magnetic site is randomly visited many times and flipped 
with a probability given by either the heat bath or metropolis algorithm in temperature steps of 
0.01 K. Both periodic boundary conditions and free boundary conditions were applied to the lat- 
tices as they were cooled and then warmed in magnetic fields of H = and 2 T. The first three 
nearest-neighbor interactions were included in the Hamiltonian with values taken from spin-wave 
dispersion measurementsQ for FeF 2 , where J x = -0.069 K, J 2 = 5.278 K, and J 3 = 0.279 K. The 
correlation functions for these three types of neighbors were calculated at each temperature step 



and averaged over ten different configurations of magnetic spins. Increasing the number of MC steps 
per temperature step dramatically did not change the results once a minimum number of steps is 
performed so that the simulation stays in quasi-static equilibrium at each temperature. The propor- 
tionality constant between C m and d(An)/dT is foundEdO to be the same in Feo.93Zno.07F2 and 
Feo.46Zno.54F2 such that, AC m — d(An)/dT, where A = 9.17 x 10~ 6 . Our sample i^eo. 25-^0.75-^2 
is near the percolation threshold concentration and cannot obtain long-range order associated with 
J2. However there is short-range order which produces a rounded peak in the C m data at about 
T = 23 K. d(An) / dT vs. T data agree with the C m data near this peak when the proportionality 
constant A — 9.17 x 10~ 6 is used as shown in Fig. 1. However, at low temperature C m and d(An)/dT 
are no longer proportional. 

MC simulations provide the opportunity to investigate the relative importance of the different 
nearest-neighbor interactions on the birefringence. The second-nearest-neighbor interaction domi- 
nates the energy Hamiltonian. Therefore the specific heat is dominated by the temperature derivative 
of the second-nearest-neighbor correlation at the percolation threshold concentration. However the 
simulations indicate that d(An)/dT is significantly governed by the temperature derivative of the 
third nearest-neighbor correlation. This is evident because the two different correlations have con- 
trasting low temperature behaviors, with only the third-nearest-neighbor one changing the sign of 
its slope at low temperature as shown in Fig. 2. 

In Fig. 3 we show the MC results for the Ising C m in the percolation threshold concentration 
simulation along with the d — 1 Ising and Heisenberg exact models. It is reasonable to expect 
that the d = 1 Ising model should correspond roughly to the behavior of the simulation. In both 
systems the spins have an average of two neighbors, there is a similar energy gap and each should 
have a vanishing specific heat at high temperature. The d = 1 Heisenberg model, on the other 
hand, has a finite specific heat at zero temperature (the model is not strictly followed by any real 
system at low temperature) . This comparison does help to explain the fact that the MC data, shown 
again in Fig. 4 as the small amplitude curve, do not correspond very well to the experimental data. 
Although the real system has a large anisotropy, it nevertheless has a Heisenberg interaction. While 
the system exhibits an energy gap and its specific heat falls to zero as the temperature is decreased, 



a significant amount of specific heat near T = 10 K is attributable to this Hcisenberg interaction. 
Since we are interested in the relative contributions to the biregringence from the frustrating third 
nearest interaction, we choose to multiply the pure Ising behavior of the MC data by a factor of 1.5, 
thereby obtaining the good fit for T < 20 K exhibited in Fig. 4. 

As mentioned above, at magnetic concentrations x = 0.93 and x = 0.46 it has been shownEJ'EJ 
that the proportionality between C m and d(An)/dT is approximately the same as for pure FeF-2. If 
we assume the same proportionality for the percolation threshold concentration in the temperature 
regime where C m and d(An)/dT are proportional and that the third nearest-neighbor correlation 
dominates the behavior of the birefringence, we obtain 

AT V r r d(< SiSj >2) -T V r r ^ (< SiS > >3) (K\ 

i<j i<j 

Furthermore, in this temperature regime we expect 

y^ ejCj < SjSj >2= 2 \] ej€j < SjSj > 3 , (6) 

i<3 i<3 

where the factor of two arises because there are eight second-nearest neighbors and only four third- 
nearest neighbors. This allows us to calculate I3 in the previous equation to be 4.85 x 10~ 5 . With 
this value of I3 and the third nearest-neighbor correlation from MC simulations we obtain a good fit 
to d(An)/dT at low temperature using the same factor of 1.5 used for C m , as shown in Fig. 5. We 
stress that we have calculated I3 in a regime where C m and d(An)/dT are proportional and used 
this value to fit the simulation results to d(An)/dT, so it is striking that the fit works well below 
T = 8 K where C m and d(An)/dT are no longer proportional. This provides strong evidence that the 
third nearest-neighbor correlation dominates the birefringence at low temperature and in the region 
where C m and d(An)/dT are proportional. However a small contribution from other correlations 
can not be ruled out. Detailed fits between the correlations calculated by MC simulation and the 
birefringence below 8 K suggest such contributions are less than ten percent. 

To understand why the small, frustrating third-nearest-neighbor interaction affects the third- 
nearest-neighbor correlation but has little effect on the second-nearest-neighbor correlation, we turn 
to the computer simulations. When the simulation is run with only the Ji interaction included in 
the Hamiltonian, all three nearest-neighbor correlations have the same temperature dependence as 



expected. Close inspection of the MC simulation result show that 1.60% of the sites flip at low 
temperature as a result of the frustration. This is close to the probability of having a site that 
contains a spin having J3 interactions but no J2 interactions which is 1.71%. This suggests that it is 
the spins that are not dominated by the J2 interactions that can order via the J3 interaction. The 
ordering of these spins is thereby dominated by the J3 interaction and the temperature derivative of 
the difference between third nearest-neighbor correlations calculated with and without the J3 inter- 
action present shows a typical Schottky peak around a temperature comparable to the interaction 
strength times average number of neighbors. Since the probability of finding such spins decreases 
dramatically above the percolation threshold concentration, no such effect should be observed at 
larger concentrations. 

We finally turn to the effect of an applied field on the behavior of d(An)/dT and, hence, the 
third-nearest-neighbor correlations. Having determined the zero-field proportionality between the 
birefringence data and simulations of the J3 correlations at low temperature, we can now test whether 
the same proportionality holds in an applied field even though the correlations themselves change 
significantly with an applied field at low temperatures. In Fig. 5 we show the measured behavior 
of d(An)/dT vs. T in an applied field H = 2 T along the c-axis. Also shown in the figure is the 
MC prediction for the behavior based on the value of I3 obtained for zero field. The fit is clearly 
good. Hence, the value of I3 has not changed significantly and the simulation correctly describes 
the temperature dependence of the birefrincence in the presence of an applied field. This result is 
important since the birefringence technique has been used extensively^ to study the RFIM which is 
realized by applying a field to dilute anisotropic antiferromagnets above the percolation threshold 
concentration. Near a phase transition the proportionality between C m and d(An)/dT holds in 
a field and the birefringence technique will yield the proper critical behavior. This is important 
because the birefringence technique invariably yields higher quality data. 

We have demonstrated that a small frustrating interaction, J3, affects the corresponding corre- 
lation, causing a reversal in the sign of the temperature derivative at temperatures comparable to 
the interaction strength J3. In turn, this causes a reversal in the temperature derivative of the 
birefringence in the case of Feo. i^Zn®. 75F2. The behavior of the birefringence in this system is 



fortuitious because it allows us to determine that third-nearest-neighbor correlation is the dominant 
one responsible for the birefringence. In most other cases of d = 3 systems, one cannot determine the 
source of the birefringence because all the correlations have similar temperature dependences. The 
fact that the proportionality, given by A, is nearly constant for x between the percolation threshold 
and pure concentrations in Fe x Zn\- x F2 strongly suggests that the birefringence is dominated by 
the third-nearest-neighbor correlations for all concentrations. We used the value of J3 measured 
for pure FeF2- The good agreement between the MC simulation and the data suggests that the 
strength of the frustrating interaction J3 is close to the same value in the percolation threshold 
concentration sample. Spin-glass-like behaviorEil near the percolation threshold has strongly hinted 
at frustration in FeQ.i^ZnQj^Fi-, but the strength of the frustration has not been previously deter- 
mined. Far above the percolation threshold concentration, or at high temperatures near percolation, 
this frustrating interaction has little affect on the ordering processes or correlations, which all have 
the same temperature dependence. 

By comparing experimental data for C m and d(An)/dT and interpreting the results using MC 
simulation results, we have been able to explain why the proportionality between C m and d(An)/dT 
fails at low temperature near the percolation threshold concentration but holds at higher temper- 
atures. From the modeling we have done, we have added to the understanding of the origins of 
the birefringence in the experimental system. We also have shown that the birefringence technique 
is well behaved even in a field and, as is well documentecD, the critical behavior of the RFIM is 
faithfully characterized with the birefringence technique. 

This work was funded by Department of Energy Grant No. DE-FG03-87ER45324. 
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FIG. 1. The magnetic specific heat C m /R = C v /R — B, where B is the background, versus T for 
Feo.25Zno.75F2- The phonon contribution to the specific heat has been subtracted as discussed in the 
text. Also shown is the temperature derivative of the birefringence d(An)/dT versus T scaled with the 
proportionality found at other concentrations. 
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FIG. 2. The first three nearest-neighbor correlation functions T; = J^< ■ > £i£j < SiSj >i /N, where iV 
is the number of sites, versus T for Feo. 2sZno. 75F2 from MC simulations. Note the change in sign of the 
slope in the H — T curve at 8 K. 
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FIG. 3. The magnetic specific heat C m /R versus T for various models. The C m per spin predictions have 
been divided by four in anticipation of a comparison with Feo.25Zno.75F2, where the molar specific heat is 
per Avogadro's number of molecules, which includes the zinc and is four times larger than the number of 
spins. 
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FIG. 4. The magnetic specific heat C m / R = C p /R — B, where B is the background, versus T for 
Feo.25Zno.75F2- The phonon contribution to the specific heat has been subtracted as discussed in the 
text. The lower thin line is a fit using only the second nearest-neighbor correlation function from MC sim- 
ulations. The MC results are then scaled by a factor of 1.5 to obtain a good fit to the experimental data at 
low temperature. 
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FIG. 5. The temperature derivative of the birefringence d(An)/dT versus T for Feo.25Zno.75F2- Note 
the change in sign in this H = T curve at 8 K. The curve through the H — data is a fit using only 
the third-nearest-neighbor correlation function from MC simulations with the coefficient I3 and then scaled 
by the same factor of 1.5 used to fit the C m data in Fig. 4. The curve through the H = 2 T data is the 
corresponding fit using the third-nearest-neighbor correlation from the H = 2 T MC simulation. 
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